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ABSTRACT 

We investigate magnetic field amplification in a turbulent velocity field with 
nonzero helicity, in the framework of the kinematic Kazantsev-Kraichnan model. 
We present the numerical solution of the model for the practically important 
case of Kolmogorov distribution of velocity fluctuations, with a large magnetic 
Reynolds number. We find that in contrast to the nonhelical case where growing 
magnetic fields are described by a few bound eigenmodes concentrated inside the 
inertial interval of the velocity field, in the helical case the number of bound 
eigenmodes considerably increases; moreover, new unbound eigenmodes appear. 
Both bound and unbound eigenmodes contribute to the large-scale magnetic 
field. This indicates a limited applicability of the conventional alpha model of a 
large-scale dynamo action, which captures only unbound modes. 
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1. Introduction 



Magnetic fields in planets and stars, protogalaxies and galaxies, and possibly intergalac- 
tic medium are generated due to random stretching of magnetic field lines by t urbulent mo- 
tion of highly conducting fluids or plasmas in which these lines are frozen fe.g..lLvnden-Bell 



1994 ; |ParkeJ 
2006 



1979 



Moffatt 



Kulsrud fc Zweibe] 



19781: lKulsrudll2005l : IZweibel fc Heilesll 19971 : ISchekochihin fc Cowlev 



20071 ). It is natural to expect that such a dynamo mechanism can 



amplify magnetic flelds at scales smaller than the correlation scales of the velocity flelds. 
However, magnetic flelds observed in astrophysical systems often appear to be correlated at 
larger scales. Such magnetic flelds can be explained if one assumes that the velocity fleld 
v(x, t) possesses nonzero kinetic helicity, that is, if = J v ■ (V x v) d^x ^ 0. To describe the 
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large-scale fields, one genera lly applies the alpha model (ISteenbeck. Krause &: Radlerlll966l : 
Moffattlll978l : lKulsrudll2005l ). This model is obtained if one averages the induction equation 



5tB = V X (v X B) + r/V^B, (1) 

where rj is resistivity or magnetic diffusivity, over small-scale fluctuations of the velocity and 
magnetic fields, assuming that these fluctuations are much weaker and concentrated at the 
scales much smaller than the scales of the growing large-scale field. As a result one obtains 
the alpha model equation for the large-scale magnetic field B(x, t), that is, the magnetic 
field averaged over the scales larger than the scales of velocity fluctuations. 



dtB = aVxB + (3V^B. 



(2) 



In this equation, a ~ v ■ (V x v)ro, and (3 ~ VqIq, where Vq is a characteristic velocity, /q is 
the characteristic scale, and tq ~ Iq/vq is a characteristic eddy turnover time (decorrelation 
time) of fluid fluctuations. In astrophysical systems resistivity is very small, so that turbulent 
velocity field effectively generates small-scale magnetic fluctuations, which grow much faster 
than the large-scale field. In this case, the applicability of the mean-field equation (121) is 
questionable (see, e.g., the disc ussion in IVainshtein fc Cattaneolll992l : lGruzinov fc Diamond 
1994J : Isiackman fc Field! 120021 ). 



To investigate this question, we use a solvable model of kinematic dynamo, introduced by 
Kazantsev (1968) and Kraichnan (1968). In this model, velocity is assumed to be Gaussian, 
with zero mean, (v) = 0, and the covariance tensor 



(x,t)t;^'(x',t')) = /t''(|x - X |)(5(t - t') 



where «;*■' is an isotropic tensor. 



k'H^) 



t^N 5^' - 



x^x^ 



x^ 



+ t^L- 



XX- 



2 +ge'^^x^. 



X 



(3) 



(4) 



Here () denotes ensemble average, e^^'^ is the unit anti-symmetric pseudo-tensor and sum- 
mation over repeated indices is assumed. The first two terms on the right-hand side of dl]) 
represent the mirror-symmetric, nonhelical part, while function g{x) describes the helical 
part of the velocity fluctuations. For an incompressible velocity field (the only case we are 
considering here), we have hn{x) = i^l{x) + xfi;'^(x)/2, where the prime denotes derivative 
with respect to x = |x|. Therefore, to describe the velocity field, we need to specify only 
two independent functions, say, kl{x) and g{x). The magnetic field correlator can similarly 
be expressed as 



{B\^,t)B\0,t)) = M^[6'^- 



x'^x^ 



x^ 



+ Mi 



'^ nf3 



XX- 
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+ Ke'^^x^, 
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where the corresponding solenoidality constraint implies Mn{x, t) = Ml{x, t) + (a;/2)M^(a;, t). 
To describe the magnetic correlator we therefore need only two functions, ML{x,t) and 
K{x,t), corresponding to magnetic energy and magnetic helicity. 

Suppose that the velocity field ([3]) is given. The problem is then to find the correlation 
funct ion ([5l) of the magnetic field. In the nonhelical case, g{x) = 0, it was established 



by iKazantsevI (119681 ) that the problem is reduced to a quantum mechanical problem with 
imaginary time: 



- dt^p = n^p. 



(6) 



More precisely, given the velocity correlator (^ and magnetic resistivity t], one constructs the 
self-adjoint Hamiltonian 7i. The magnetic correlator ML{x,t) is then mapped to the "wave 
function" ip^x, t). If the equation IQ has growing solutions ipi^x, t), corresponding to negative 
eigenvalues of 7i, then dynamo action is possible. One can show that the smaller the resistiv- 
ity f], the deeper the effective potential in the quantum mechanical problem ([6]), and, there- 
fore, nonhelical d ynamo is always possible when t he magnetic Reynolds number, Rm oc rj~ ^ 



is lar ge enough (IVainshtein fc Kichatinovl Il986l : iBoldyrev fc Cattaned l2004l : llskakov et al. 



20071 ). The quantum mechanical representation ([6]) is important since it ensures that the 
eigenvalues are real and the eigenfunctions are mutually orthogonal, which allows one to ap- 
ply a variational principle for estimating dynamo growth rates. The Kazantsev-Kraichnan 
model of nonhelical dynamo action has been well investigated in the literature. 

The situation is more complex when the velocity field possesses nonzero helicity, i.e., 
g{x) j^ 0. In this case, given kinetic energy k.l{x) and kinetic helicity g{x), one needs 
to solve two coupled partial differential equations for functions Mi{x^t) and K{x,t) re- 
lated to magnetic energy and magnetic helicity. Such equations were first derived by 
Vainshtein fc Kichatinovl (119861 ) . Due to their complexity, there have been relatively fewer re- 



sults obtained for the helical case (e.g., 



1992 



Kim fc Hugheslll997l : lBrandenburgl 



Meneguzzi. Frisch fc Pouauet 



2001 



Blackman fc Field 



1981 



2002 



Kulsrud fc Anderson 



Brandenburg fc Subramanian 



20051 ). However, it is the helical case that is practically more important since astrophysi- 
cal systems generally possess nonzero helicity. Moreover, while magnetic fields at small 
(velocity) scales are naturally expected on the base of equation ([6]), explanation of astro- 
physically observed large-scale magnetic fields commonly requires helicity effects, as, e.g., in 
equation ([2]). 



Recently, it has been established in lBoldyrev. Cattaneo fc Rosnen (120051 ) that the I Vainshtein fc Kichati 



(119861 ) equations also possess a self-adjoint structure, and can be reduced to a quantum me- 



chanical "spinor" form: 



- (9i^" = n'^''^p 



apjfi 



(7) 
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where a = {1, 2} and summation over repeated indices is assumed. Similarly to equation ([6]), 
the self-adjoint Hamiltonian 7i°^ depends on kinetic energy and kinetic helicity, kl{x) and 
g{x), and on magnetic resistivity r/. The two components of the il)°'{x,t) function are then 
related to functions Mi^x^t) and K{x,t) in the magnetic correlator (^. 

Similar structures of the dynamo equations ([6]) and ([7]) allow one to investigate them 
on the same footing. In particular, one can compare the growth rates and the eigenvalues 
associated with helical and nonhelical dynamo action, and address the important question 
of whether the large-scale magnetic field generated due to helical dynamo action is described 
by the alpha model ([2]). This is the goal of the present work. We assume that the ve- 
locity field has the Kolmogorov spectrum and that the Reynolds number and the magnetic 
Reynolds number are large. Then we present the full numerical solution of the helical dynamo 
model d?]), i.e., we find its eigenvalues and eigenfunctions. We demonstrate that at least at 
the kinematic stage, the alpha model ([2]) may provide a nonadequate description of large- 
scale fields, since it misses the rapidly growing large-scale eigenmodes. In the next section 
we present our main results; a detailed discussion will be presented elsewhere, nonhelical 



2. Numerical solution of the helical dynamo model 



The self-adjoint equations for functions Ml{x, t) and K{x, t) were derived in lBoldyrev. Cattaneo fc Ros 



(120051 ). We will be interested in eigenmodes of these equations, and therefore assume that 
both functions depend on time as exp(At). It is then convenient to introduce the auxiliary 
functions W2{x) and w^i^x) defined as: 



Mr 



V2e 



xt 



^xt 



x^ 



-W2{X) 
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[x^U73(x)]' 



The eigenmode equation that we solve then takes the form: 



X X 

dx x^ 



\f2fi d 2 



dx 
d_ 
dx x'^ dx 



^2_d__B _d_^2 _ ^ 



W2 
W3 



(9) 



where 



E 

A{x) 
b\x) 
C{x) 



2 dr dr 

Ai LhJu/ \JbJu 



v^[2r/ + fi:^(0)-fi:^(x)], 
2?7 + Ki:(0) -/tL(x), 
^f2\g{^) - g{x)\x, 



(10) 



and primes denote derivatives with respect to x. Equations (l8|)- flT0|) describe the growth of 
the magnetic field in the Kazantsev-Kraichnan model. Equation ([9]) is self-adjoint, which 
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guarantees that all growth rates A are real. The system (Q cannot be solved analytically 
for general velocity correlation functions kl{x) and g{x). Below we solve equation iQ nu- 
merically, and concentrate on functions W2{x) and 1^3 (x) since they define the magnetic field 
correlator uniquely and contain all the information about magnetic energy and magnetic 
helicity. 



It is useful to note the Fourier transformed version of Equation 



K'^{k) 



F{k) 1 6'^ 



A;2 



+ iG{ky^^k\ 



(11) 



where functions F{k) and G{ k) can be expressed in t erms of the three-dimensional Fourier 
transforms of kl{x) and g{x) (JMonin fc Yaglomlll971l ). and 



{B\\^,t)B*^%t)) = Fsik.t) (5^^ 



k'y 






(12) 



where Fsi^k, t) is the magnetic energy spectral function, (|B(k, t)p) = 2FB{k, t), and Hsik, t) 
is the spectral function of the electric current helicity, (5* (k, t) ie^^^y B^{k,t)) = HB{k,t). 

First, as easily derived from ([9]) the asymptotic behavior of functions W2{x) and Ws{x) 
when X -C ■\/?7/k'^(0) is 



W2 = x'^ + 0{x^ 



ws = ^x^ + 0{x^ 



(13) 



where, without loss of generality, we use scahng W2/x'^ — > 1 as x ^ 0, and coefficient C, is 
a free parameter, related to the averaged helicity of the electric current, (B ■ (V x B)) = 
6K{0) = -15^2^. 

Second, the asymptotic behavior of functions W2{x) and ^3(0:) as x — > 00 depends on 
the magnetic field growth rate A. If A > Aq = g^{0)/[KL{0) + 2r]], then as x — >• 00 the 
asymptotic eigenfunctions are 



W2 
W3 



oc e 






cos(fcjX + 0) 
sin(A;jX + 0) 



(14) 



where kr = y/X — Ao/A//ti(0) + 2r/ and ki = \^\)/ ^/f^iiO) + 21]. For A < Aq, the asymptotic 
behavior of eigenfunctions W2{x) and W3(x) at x ^ 00 becomes a mixture of cosine and 
sine standing waves, cos{kfx + 0^) and sm{kfx + 0^), with wavenumbers kf = {y/Xo ± 
a/Ao — A)/a/kl(0) + 2?7. This asymptotic behavior suggests an analogy between equation ([9]) 
and quantum mechanics. The spatially localized (bound) eigenfunctions (TT^ with A > Aq 
correspond to "particles" trapped by the potential provided by velocity fluctuations, while 
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the nonlocalized (unbound) eigenfunctions with A < Aq correspond to "traveling particles." 
Eigenvalue A = Aq corresponds to the fastest growing unbound eigenmode. 

To find the eigenmodes, we proceed as follows. For given values of A and ^ we integrate 
ordinary differential equations iQ numerically by the fourth-order Runge-Kutta method. 
Due to quite disparate scales present in the problem, we use a nonuniform numerical grid 
with the grid steps chosen as Ax{x) oc B^^'^{x) for x > y/rj and Ax{x) oc {x/ ^)B^f'^{^) 
for X < y/rj. The integration is done starting at a small value of x = Xmin ^ \/vT^^iW)i 
where the asymptotic result (iT3l) holds. We integrate up to a large value of x = Xmax that 
is chosen in such way that the numerical solution is still stable, while it has already reached 
its asymptotic behavior for large values of xo By matching the numerically calculated 
values if2(xmax) and w^i^Xraax) to their analytic asymptotic solutions at x = Xmax, we find 
the eigenvalues (the magnetic field growth rates) A and coefficients ^. It turns out that 
these eigenvalues are discrete for the bound (localized) eigenmodes, i.e., A = A„ > Aq and 
i = in in = 1, 2, 3...). In contrast, the eigenvalues of the unbound (nonlocalized) modes are 
continuous; they exist for any choice of A < Aq and for any realizable value of ^o 

To study a realistic case, we consider velocity correlation tensor (llip with the Kol- 
mogorov power velocity spectrum. It is important to note that in the Kazantsev-Kraichnan 
model ([3]), the velocity field enters the eigenfunction equations ([6], [71 [9]) only in the form of 
turbulent diffusivity K^-'dx — x'|) = J(f*(x, t)f-'(x', t'))(i(t — t'). I n the Kolmo gorov turbu- 



lence, the latter scales as vil ~ Z^^'^, where / = |x — x'| (see, e.g., lFrischlll995l ). Comparing 
this with formula ([3]), we find that the Kolmogorov scaling implies k(x) ?a k(0)(x//o)'^^'^ ~ 
voIq{x/IqY/'^. Without loss of generality we take Iq ~ 1, fo ~ 1, and therefore 

F{k) = k-^^'^, 

G{k) = -hk-^Fik), 2<A;<A;^ax. (15) 

Here the minimal cutoff wavenumber kj^in = 2, the maximal cutoff wavenumber /cmax ~ 
2[kl(0)/u]^^^ ^ 4z/~^/^ is determined by the plasma kinematic viscosity z/, and the helic- 
ity parameter h must satisfy the realizability condition —l<h<l; the velocity field is 
maximally helical when \h\ = 1. The resulting growth rates A„ of the bound (localized) 
eigenmodes. An > Aq, are shown in Figured! where the magnetic diffusivity is chosen to be 
7] = 10~^. The growth rates are mesured in units of large-scale eddy turnover rate ~ Vq/Iq. 



^ Our results do not depend on the exact choice of the boundaries of the computational interval a::min ^ 

-^ _ ■^max- 

^ Given Mi^(x,t), function K[x,t) cannot be chosen arbitrarily, its Fourier image must satisfy the realiz- 
ability condition \Hf!i k,t)\ < Ff!(k, t). Analogously, given kl{x), function g{x) must also satisfy a similar 
realizability condition ( Moffattill978[ ). The later results in the condition — 1 < /i < 1 in Eq. ([TS]) . 
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The linear-scale plots A and B correspond to h = 1 and h = 0.1 respectively, while fcmax = 3, 
which is consistent with the Reynolds number being of order unity (i.e., a single-scale ve- 
locity field). The logarithmic-scale plots C and D correspond to h = 1 and h = 0.1, while 
^max = 3000, which is consistent with large Reynolds and large magnetic Prandtl num- 
bers. Finally, the logarithmic-scale plots E and F correspond to h = 1 and h = 0.1, while 
kmax = 3 X 10^, which is consistent with a very large Reynolds number and a small Prandtl 
number. 



3. Discussion and conclusion 

We find that when the Reynolds number is of order unity (a single-scale velocity field), 
the magnitude of the kinetic helicity parameter h does not have much effect on the bound 
magnetic eigenmodes (plots A and B in Fig. [1]). However, in the case when the Reynolds 
number is large and the velocity fluctuations extend over a large range of scales, the bound 
eigenmodes are signiflcantly affected by kinetic helicity at scales larger than the viscous 
scale. (Magnetic fluctiations at the scales much smaller than the viscous scale, which are 
excited in the case Pr >> 1, are not signifi cantly affected by magnetic helicity, which is 
consistent with previous considerations; e.g.. iKulsrud fc AndersonI (119921 ).) The number of 



bound eigenmodes increases considerably when the kinetic helicity increases, and the cor- 
responding eigenvalues, A„, become strongly concentrated near the eigenvalue of the fastest 
unbound eigenmode, Aq. This result follows from a nearly uniform distribution of A„ on the 
logarithmic-scale plots in Fig. [1] 

Moreover, we observe an important fact that in all the cases in Fig. [T] the growth 
rate of the first bound eigenmode, Ai, happens to be very close to the growth rate of the 
fastest growing unbound eigenmode, AqIj We conjecture that for high Reynolds numbers 
the potential in Qj always has a shallow bound state Ai (such that Ai — Aq <C Aq) whose 
characteristic scale is much larger than the velocity correlation scale, due to flT^ . This has 
an important consequence for the dynamo mechanism. To understand it, we note that the 
alpha dynamo model (jSj), describing the large-scale magnetic field, does not capture the 
bound eigenmodes. However, our analysis shows that shallow bound eigenmodes have faster 
growth rates and large correlation lengths, so at a given scale x they may rapidly become 
dominant over the unbound modes. In practical applications, this means that such modes 
rather than the modes described by ([2]) become essential in the large-scale magnetic field 



■^We found same result for all other high-Reynolds cases that we investigated (not reported here), with 
different helicity parameters. 



configurations. In this case the conventional alpha-dynamo model ([2]) leads to an inadequate 
description of the large-scale magnetic field even at the kinematic stage of dynamo action. 
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Fig. 1. — Growth rates A„ > Aq of the bound magnetic eigenmodes for the Kolmogorov 
velocity spectrum given by Eq. (ITSl) and r/ = 10~^. Plots A and B are for h = \ and 0.1, 
while fcmax = 3 (Reynolds number Re ~ 1). Plots C and D are for h = 1 and 0.1, while 
^max = 3000 (Re ^ 1 and Prandtl number Pr ^ 1). Plots E and F are for h = 1 and 0.1, 
while fcmax = 3 X lO'^ (Re > 1 and Pr -C 1). 



